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Abstract 

Using the Lie's infinitesimal method we establish that the Dirac equation in one 
variable is integrable by quadratures if the potential V(x) is a solution of one of the 
equations of the stationary mKdV hierarchy. 

Consider the eigenvalue problem for the Dirac operator C = io~i-^ — V(x)<T2 

du -* 
(£-\)u = ia 1 —-(V(x)a 2 + \)u = (1) 
dx 

where o"i,cr 2 are 2x2 Pauli matrices, u— (u\(x), u 2 (x)) T , V(x) is a real-valued function 
and A is a real parameter. We remind that ([[]) is one of two equations composing the Lax 
pair for the mKdV equation 

v t + v xxx - 6v 2 v x = (2) 
integrable by the inverse scattering method (see, e.g. [I], 0). Next, as the identity 

holds, components of the vector-function u fulfill the stationary Schrodinger equation 

'' 2 " l + ((-l) i+1 ^-V 2 +Au i ^0, i = l,2. (3) 



dx 2 \ dx 
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The aim of the present paper is to show that there exists an initimate connection 
between integrability of system (JJ) (in what follows we will call it the Dirac equation) by 
quadratures and solutions of the stationary mKdV hierarchy. 

Integrability of system (JJ) will be studied with the use of its Lie symmetries. As usual, 
we call a first-order differential operator 

X = S(x)-fj- + ri(x), 
ax 

where £ is a real-valued function and rj is a 2 x 2 matrix complex-valued function, a Lie 
symmetry of system ([!]) if commutation relation 

[£, X] = R(x)C (4) 

holds with some 2x2 matrix function R(x) (for details, see, e.g. j3|). 

A simple computation shows that if X is a Lie symmetry of system (H), then an 
operator X + r(x)C with a smooth function r(x) is its Lie symmetry as well. Hence we 
conclude that without loss of generality we can look for Lie symmetries within the class 
of matrix operators X = r)(x). Furthermore, inserting X = rj(x) into (f|) and computing 
the commutator yield that the matrix rj(x) is necessarily of the form 



T) 



f(x) g(x) 
h{x) -f(x) 



(5) 



where f(x),g(x), h(x) are arbitrary solutions of the following system of ordinary differen- 
tial equations: 

4L = i\( g -h), ^ = 2t\f + 2gV, ^ = -2i\f - 2hV. (6) 

With a solution of system @ in hand we can integrate the initial equations ([!]) by 
quadratures using the classical results by Elie Cartan [[|. Since these results are well- 
known we will give them without derivation in the form of the following lemma. 

Lemma 1 Let the functions f(x),g(x), h(x) satisfy system ^). Then the general solution 
of the Dirac equation is given by the formulae 

Ul (x) = C l {R{x) + f{x)){h{x))-^{R\x)-X)- l l\ 

u 2 (x) = Ci (h(x)) 1/2 (R 2 (x) - A)~ 1/2 , 1 ' 
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where A = f 2 {x) + g(x)h(x) is constant on the solution variety of system (|6[), 

dx \ 



R{x) 



VA tanh \C 2 - iX^/A J 
y^A tan ^C 2 + iXV^A J 



dx 

g(x) i 



under A > 0, 
under A = 0, 
under A < 



and Cx,C 2 are arbitrary complex constants. 



However, solving system of ordinary differential equations ((^) is by no means easier 
than solving the initial Dirac equation. This is a common problem in applying Lie sym- 
metries to integration of ordinary differential equations. The key idea of our approach is 
to restrict a priori the class within which Lie symmetries are looked for and suppose that 
they are polynomials in A with variable matrix coefficients. 

Introducing the new dependent variables tjji(x) , t/j 2 (x) 



g(x) 
h(x) 



dip\ 



+ 2Vij 2 



i 

4A 
I 
2 

~(ipi(x) -ip 2 ( x )) 



dx 

(ipi(x) + ijj 2 (x)), 



we rewrite (0) in the following equivalent form: 



d 2 i) x 
dx 2 



dip 2 



dibo dV 

-AX 3 ^ 1 + 2V-p- + 2^ 2 — , 

dx dx dx 



2Vifj 1 



(9) 



As mentioned above solutions of system (^j are looked for as polynomials in A, namely 

n n 

<M*) = £^)(2A) 2fc , V^) = E^)(2A) 2/C . (10) 
fe=i fc=i 

Inserting the expressions ([It]) into (|9[) and equating the coefficients by the powers of 
A yield p n = and 

dr k 
dx 



2Vp kl k = 1, 



,n 



d 2 Pk dr k dV 

-TT = 2V— + 2— r k -pk-u k = l,...,n 
dx z dx dx 



'12) 



where we set by definition = 0. Eliminating from (|TTD, (JT^) the functions r^x) we 

get recurrent relations for the functions Pk{x): 

{ d 2 dV 1 

Pk-i(x) = + ±-faD- x V + AV 2 Pk (x), k = n, n - 1, . . . , 0. (13) 



Q 

Here D" 1 denotes integration by x. 

A reader familiar with the theory of solitons will immediately recognize the operator 
Q as the recursion operator for the mKdV equation (|2|) (see, e.g. @, §)• Acting re- 
peatedly with this operator on the trivial symmetry So = yields an infinite number of 
higher symmetries Si, S 2 , ■ ■ ■ admitted by the mKdV equation j5|. Hence it is not difficult 
to derive that the functions Pk,k = 0, . . . ,n — 1 are linear combinations of the higher 
symmetries S±, . . . ,S n with arbitrary constant coefficients C{ 

k 

Pn-k{x) = X) CA+i-i, k = l,...,n, (14) 

i=l 

where Si are determined by the recurrent relations 

Sl+l(x) = d^T + A lx'i V ^ S ^ dy + AV2S ^ x )i < = 1, ■ ■ ■ , n - 1 

x 

with S, d ^ f f . 

The above formulae (13) give the general solution of the first n equations from fli"3|). 
Inserting these into the last equation yields equation for the function V(x) of the form 

n+l 

CkSn+2-k = 0. (15) 

k=l 

As 5*i = (|1|) is nothing else than an equation of the stationary mKdV hierarchy, which 
is obtained from the higher mKdV equations by setting v(t, x) = v(x + Ct), C = const. 
Integrating equations fllTf ) yields 

k x 

r k (x) = 2j2Q J V(y)S k+1 ^(y)dy + C k , k = l,...,n, (16) 
where Ci are arbitrary complex constants. 
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Thus, the formulae fllCf), fll4|), (p)|), ( |T6D give the general solution of the system of 
determining equations (|TT|) , (|l"2D within the class of functions of the form fllPl) . This 
means, in particular, that provided the function V(x) is a solution of equation ( |H)D with 
some fixed n and C\, . . . ,C n , the Dirac equation possesses a Lie symmetry. Hence we 
conclude that it is integrable by quadratures due to Lemma 1. Consequently, we have 
proved the following remarkable fact. 

Theorem 1 Let V(x) be a solution of an equation of the mKdV hierarchy of the form 
(175[). Then the Dirac equation (|7[) is integrable by quadratures. 



Note that the equations of the stationary mKdV hierarchy are transformed to the 
equations of the stationary KdV hierarchy with the help of the Miura transformation and 
the latter are integrated in ^-functions with any n G N 0. 

In conclusion let us demonstrate how the above procedure works for the simplest case 
n = 1 . With this choice of n equation ( |15D reads 

C 2 dV _d 3 V QV 2 dv () (17) 
C\ dx dx 3 dx 

which is exactly the stationary mKdV equation and is obtained from (§) via the Ansatz 
v (t,x) = V(C 2 x-C\t). 

A simple computation yield the form of the coefficients of the Lie symmetry (^]) 

d 2 V \ 
^-2C 1 V*-C 2 -4C 1 \ 2 )), 

_C l V 2 -^C 2 ~2C 1 X^j , (18) 

+ dV 2 + X -C 2 + 2dA 2 ^ 

which satisfy the determining equations (|6|) inasmuch as the function V(x) is a solution 
of the stationary mKdV equation. 

Thus, the Dirac equation with potential V (x) satisfying the stationary mKdV equation 
(|17D is integrable by quadratures and its general solution is given by formulae (0), 



fix) 






g(x) 


1 

2 


( 1 dx 


h(x) 


1 
2 


I 1 dx 



Note that due to the remark made at the very beginning of the paper components of 
the function u fulfill the stationary Schrodinger equation (|3]). This is in a good accordance 
with results of papers || ||. 
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